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ABSTRACT 


Light curves of microlensing events involving stellar binaries and planetary systems 
can provide information about the orbital elements of the system due to orbital modu¬ 
lations of the caustic structure. Accurately measuring the orbit in either the stellar or 
planetary case requires detailed modeling of subtle deviations in the light curve. At the 
same time, the natural, Cartesian parameterization of a microlensing binary is partially 
degenerate with the microlens parallax. Hence, it is desirable to perform independent 
tests of the predictions of microlens orbit models using radial velocity time series of the 
lens binary system. To this end, we prese nt 3.5 years of RV mon itoring of the binary lens 
system OGLE-2009-BLG-020L, for which Skowron et al. ( 2011 1 constrained all internal 
parameters of the 200-700 day orbit. Our RV measurements reveal an orbit that is 
consistent with the predictions of the microlens light curve analysis, thereby providing 
the first confirmation of orbital elements inferred from microlensing events. 


Subject headings: gravitational lensing: micro - techniques: radial velocities - (stars:) 
binaries: general - (stars:) binaries: spectroscopic 


1. Introduction 

While the general-relativistic microlensing effect has been repeatedly observed, very few direct 
tests of the microlensing model solutions have been possible. This is because microlensing is 
inherently a rare transient phenomenon, and the lenses observed are often extremely faint. Because 
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a microlensing event requires that two stars at different distances align with each other along the 
line-of-sight to better than ~ 1 mas, in the densest parts of the sky only about 1 in a million stars 
is exp ected to be undergoing a microlensing event at any given time (jPaczvnskilll99ll : lHan fc Gould 
19951 1. These events are relatively brief (~ 2 months) and (effectively) unrepeatable. In addition, 
since the only light required to study the gravitational potential of the foreground lens is provided 
by the background source, the lens itself can, in principle, be completely non-luminous. Typically, 
lens stars are M dwarfs that reside more than halfway to the center of the Galaxy at > 4 kpc, and 
are thus very faint, making them extremely difficult or impossible to followup after the event is 
over. This means that while the ensemble of microlensing detections is robust, very few individual 
lens stars can be studied in detail aside from what can be learned from the lensing event. More 
importantly, there have been relatively few conhrmations of the complex microlensing modeling 
process even as the number of parameters expands to include more effects. 

Most previous tests of microlensing models have focused on conhrming that the measured 


bright ness of the lens star is consistent with the predicted lens mass in microlensing. iBennett et al 


([20101) made an independent confirmation of a microlensing model based on adaptive optics obser¬ 
vations of OGLE-2006-BLG-109. They demonstrated that the observed lens flux is consistent with 
the predicted lens mass and distance made by measuri ng the parallax and the hnite source effects 
in the light curve. Likewise iPietrukowicz et al.l (j2005l i discovered a transient in M22 and classi¬ 
fied it as likely microlensing event caused by the lens belonging M22 magnifying the background 
bulge star. They gave mass estimation based on the measured Einstein time-scale, the known 
distances and glob u lar cl uster proper motion (M = O.II^qq^Mq), which was later conhrmed by 
Pietrukowicz et al.l (l2012l ) using VLT /NAGO to sh ow that the measured lens flux corresponds to 
a mass of M = O.ISMq. In addition, iGouldl ([20041 ) conhrmed the microlens parallax measured for 
MAGHO-LMC-5 by showing that t he direction o f the relative proper motion measured from the 
separation of the source and lens by lAlcock et al.l (j200ll ) was consistent with the expectation from 
the parallax. 

In this paper, we present the hrst direct test of a microlensing detection of orbital motion. 
While the orbital period at a typical Einstein ring radius of a few AU (where microlensing sensitivity 
to companions is maximized) is generally much longer than the typical lensing timescale of tE ~ 20 
days, in some cases it is still possible to observe K eplerian motioii of a binary lens system. The 
theory of orbiting binary lenses was hrst explored bv iDominikI (jl998l ). In practice, the detectability 
of orbital motion depends on having features in the light curve that are well-separated in time 
(~ 20 days) and a relatively short orbital period (< hundreds of days). The orbital motion leads 
to changes in the shape of the caustic structure, and hence the magnihcation pattern due to the 
changing separation between the components of the lens and can also rotate the caustic on the 
plane of the sky. 


MAGHO 97-BLG-41 ([Alcock et al.ll200fll ) was the hrst binary lens to show strong deviations 
from the assumption of a static binary, but th e initial interpretatio n was that the deviations were 
due to a third body (a circumbinary planet; [Bennett et al.lll999l L and it was not until a later 
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analy sis of an independent data set that the orbi ting binary soluti on was found (jAlbrow et al 
2000) • This controversy was definitively resolved by lJung et alJ (|2013l i , who conducted a combined 
fit of the data and found that direct comparison of the orbiting binary and circumbinary planet 
models clearly preferred the binary. Orbital motion is also found in systems in which the companion 
is a planet, and was in fact observed in the second planeta ry system discovered by microlensing: 
OGLE-2005-BLG-071 ( Udalski et al. 2005; Dong et al. 2009). Hence, experience has demonstrated 
that it is important to take these effects into account when htting microlensing light curves. 

At the same time, a test of an orbital motion model would greatly increase our confidence in 
the derived orbital parameters for both stellar and planetary microlenses. Introducing the orbital 
motion effect into the microlensing models by definition increases the number of free parameters, 
raising the concern that any improvement in the £t is due primarily to htting systematics or 
correlations in the data. In addition, the orbital motion parameters are known to be correlated 
with other microlensing parameters an d effects such as th e orbital parallax effect (the effect of the 
Earth’s motion on the light curve; see iBatista et al.ll201ll ) and xallarap (motion of the source due 
to a binary). Because of the correlation with parallax, a conhrmation of the orbital motion solution 
will also translate into increased conhdence in the measured mass from the parallax effect. 


While tests of the orbital motion solutions for stellar lenses with planetary companions will 
remain difficult even with 30m class telescopes, it is possible that this test could be done for binary 
star lenses because the orbital motion signal is so much larger (e.g. radial velocities O km s“^ for 
binaries rather than O m s“^ for planets). 


One system seen t o exhibit orbit a l mot ion is the binary star lens in OGLE-2009-BLG-020, 
which was analyzed by ISkowron et al.l ([201l|). They found that it was impossible to hnd a satis¬ 
factory ht to the microlensing light curve without allowing for orbital motion of the lens. From 
this analysis, they were able to place broad constraints on the Keplerian parameters of the orbit 
(reproduced in the upper, right-hand panels of Figured]). In brief, the posteriors indicate a O.84M0 
primary with an M-dwarf companion in a 200-700 day orbit with some amount of eccentricity. 


What makes this system unusual is that it is exceptionally close for a microlensing lens system, 
Dl ~ 1 kpc, such that the lens is clearly visible in the blended light. In fact, with an I magnitude 
of 15.6, the lens is brighter than the unmagnihed source. Because the system is so bright and the 
expected radial velocity (RV) signal from the lens is so large (~ 10 km s“^), it is possible to conhrm 
the orbital motion of the lens system measured from the microlensing light curve with followup 
observations. 


In this paper, we present radial velocity followup observations of the lens system and con¬ 
hrm a microlensing orbit solution for the hrst time. We begin by comparing and contrasting the 
direct observables of binary stars as seen with radial velocity and microlensing (Section [2|). The 
microlensing and RV observations of OGLE-2009-BLG-020 are presented in Section [3| A detailed 
discussion of the RV data is given in Section Hi including a novel method to use the source star as a 
wavelength reference (Section 14. 2 j) and the hnal RV solution to the orbit (Section 14.4j) . In Section 
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5l we show that this RV solution is consistent with the constraints on the orbit from Skowron et al, 


([20111), and in Section [6| we perform a joint fit to both the RV and microlensing data to find the 
final parameters of the binary system. Our conclusions are given in Section [71 


2. Orbit Kinematics: RV vs. Microlensing Observables 

Fourteen parameters are required to completely characterize the kinematics of a binary star 
orbit. These could be the six phase-space coordinates of each body at a given time plus its mass, 
or any non-degenerate set of combinations of these quantities, for example those given in Tabled] 
Single-line spectroscopic (RV) observations yield measurements nine of these 14 parameters, while 
microlensing observations can potentially yield 13. However, the natural parameterizations of these 
two characterizations are very different. Hence, before comparing microlensing predictions with RV 
measurements, it is essential to understand each parameterization. 

2.1. Radial Velocity Parameters 

RV observations are normally thought of as yielding five (out of eight) parameters that char¬ 
acterize the internal motion of the binary, plus a spectroscopic estimate of the mass of the primary. 
These are the period P, the eccentricity e, the argument of periastron cj, the time of periastron 
tperi, and the mass function Ai = (m 2 sini)^/(mi -|- m 2 )^. These are the same as the RV observ¬ 
ables except that the invariant A4 is replaced by the observable K, the RV semi-amplitude. Of 
the remaining parameters that are not measured, the longitude of ascending node (Q), is rarely of 
physical interest and is therefore usually ignored. Thus, there are only two parameters of interest 
that are not measured, the inclination i and the mass ratio q = m 2 /mi. Note that for double-lined 
binaries, q is measured, while for extreme mass ratios (e.g., planets), Ai — {qsmi)^mi. 

For RV, three center-of-mass parameters are also known, namely the measured system radial 
velocity, Vz, and the two-dimensional position on the sky (0, i.e. the coordinates of the system in 
right ascension and declination). The last two are so intrinsic to the process of measurement that 
they are not normally even considered as “measurements”. 

Table [2] summarizes the binary parameters known from RV observations. 


2.2. Microlensing Parameters 


Since microlensing can in principle measure 13 parameters, the simplest way to characterize 
these is to specify the parameter that cannot be measured: the system radial veloci ty. In addition 


micro lensing cannot distinguish between (O, cuperi) and (180 deg — fl, Wperi — 180 deg ([Skowron et al, 


201 ih . That is, the parameters that can be derived from microlensing are identical to those from 


astrometric measurements for similar reasons: namely that microlensing effects derive from the 
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Table 1. Standard, Keplerian Parameterization of a Binary 


Variable 

Units 

Meaning 

Binary Orbit Parameters: 

^peri 

days 

Time of periastron 

a 

AU 

Semi-major axis 

e 


Eccentricity 

^node 

deg 

Longitude of ascending node 

i 

deg 

Inclination 

^peri 

deg 

Argument of periastron 

mi 

Me 

Mass of primary 

m2 

Me 

Mass of secondary 

Phase-Space Parameters: 

Dl 

kpc 

Distance to system 

e 

deg 

Coordinates of the system (RA/DEC) 


mas yr“^ 

Proper motion of the system 

Vz 

km s“^ 

Radial velocity of the system 


Table 2. RV Parameterization of a Binary 


Variable 

Units 

Meaning 

RV Orbit Invariants: 



^peri 

days 

Time of periastron 

p 

days 

Period 

e 


Eccentricity 

^peri 

deg 

Argument of periastron 

Ml = {m 2 sin i)^/(mi + m 2 )‘^ 


Mass function 

(mi 

Me 

Spectroscopic mass of primary) 

Other known/measured parameters: 


e 

deg 

Coordinates of the system (RA/DEC) 

Vz 

km s“^ 

Radial velocity of the system 
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Table 3. Microlensing Parameters 


Variable 

Units 

1 Meaning 

13 Parameters of a 

Microlensing Model; 

Uq 

Oe 

Closest approach between the source and lens 

to 

days 

Time when u{t) = uq 

tE 

days 

Einstein crossing time 

P 

Oe 

Normalized source size 

tte 


Microlens parallax vector 

s 

Oe 

Projected separation of the lens components 

q 


Mass ratio between the lens components 

a 

rad 

Angle between the binary axis and the source trajectory 

7 


Normalized, projected relative velocity of the binary 


Oe 

Relative position of the lens companion along the line of sight 

Iz 

Oe 

Relative velocity of the binary along the line of sight 

Additional Known Parameters: 

e 

deg 

Angular coordinates of the microlensing event (RA/DEC) 

e. 

mas 

Angular size of the source 

TTs 

mas 

Source parallax 
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Table 4. Parameters Derived from Microlensing Parameters 


Variable 

Definition 

Units 

Meaning 

Intermediate Parameters: 



Oe 

^*Ip 

mas 

Angular size of the Einstein ring 

I^geo 

Qe/'^e 

mas yr“^ 

Geocentric relative proper motion between the source 




and lens (magnitude) 

/^geo 

TTe 


Direction of geocentric relative proper motion 

s 

(s cos a, s sin a) 

Oe 

Projected binary separation vector 

AO 

sOe 

mas 

Angular binary separation 

Mtot 

OE/inTTE) 

Mo 

Total mass of binary 

Binary Parameters: 



mi 

Mtotl (1 + ^ 

Mo 

Mass of primary primary 

m2 

mi/q 

Mo 

Mass of secondary 

Dl 

AU/ {OeT^E + TTs) 

kpc 

Distance to the binary system 

r± 

DlO 

AU 

Physical position of the binary system 

v± 

^Li-^ geo 

km s“^ 

Projected velocity of the binary system 

Ar_L 

DlAO 

AU 

Physical projected separation of the secondary 

Arz 

DlSzOe 

AU 

Relative position of the secondary along the line of sight^ 

Av_i_ 

DlsOeI 

km s“^ 

Relative, projected velocity of the secondary 

Avz 

DlIzGe 

km s“^ 

Relative velocity of the secondary along the line of sight'^ 

{Vz 



Unknown systemic radial velocity) 


“^Note that in microlensing there is a perfect degeneracy between solutions into and out of the 
plane of the sky, i.e. (Ar^, Avz) —>■ (—Ar^, —Avz)- 







time evolution of the projected positions of the two components on the plane of the sky. 


However, from the standpoint of understanding the information content of microlens binary 
solutions, the above description is a bit too simple. First, several of the microlensing parameters 
are quite unfamiliar combinations of phase-space and masses. More important, the precision with 
which these parameters can be measnred is highly variable, with some parameters measured to 
fractions of a percent and othe r s usu ally measured only to within a factor of a few (with notable 


exceptions in I Shin et al.l l 201 ll . [2013) • Properly understanding how RV and microlensing can be 


compared requires taking these differences into account. 


First we consider the parameters that are required to characterize a caustic-crossing binary 
microlensing event, which can then be transformed into 13 physical parameters of the binary. Seven 
parameters describe the basic microlensing event. Three of these describe the nnderlying point-lens 
event, i.e., the time of maximnm to, the impact parameter uq (in units of the angular Einstein radius 
0 e)) and the Einstein timescale ^e- There are three basic binary lens parameters, the mass ratio q 
and the vector projected separation s of the companion (in nnits of 0 e) relative to the direction of 
lens-sonrce relative proper motion in the geocentric frame /.tgeo- This vector is frequently expressed 
as s = (s cos a, s sin a), where a is the angle between the binary axis and the sonrce trajectory. 
Einally, there is the normalized source size p = which is required to described the extended 

duration of the caustic crossing due to the finite angnlar size of the source 0*. Of these 7 parameters 
{to,uo,tE,q-s, p) all but two {to,uo) enter the 13-parameter binary characterization. 


The accelerated moti on of Earth can induce sufficient distortions in the lightcnrve to measure 
the “microlens parallax” (|Gouldlll992l l. 


_ ^rel Rgeo 

TTE = -7-: 

Pgeo 


( 1 ) 


where vTi-ei = AU(I1^^ — is the lens-source relative parallax (see Gould &: Horne ( 2013 1 for a 
didactic explanation). 


Becanse microlensing observations are short compared to an orbital period, the orbital param¬ 
eters are naturally formulated in terms of Gartesian phase-space coordinates, rather than Kepler 
invariants (as for RV). In addition, because the microlensing event is entirely governed by the pro¬ 
jected motion of the binary, the most robustly measured parameter of the motion is the projected, 
relative velocity of the binary, Av_l. This is parameterized as instantaneous rates of change of s 
and a respectively, which yield 7 = (7||,7_l) = {ds/dt/s,da/dt). 


The remaining two parameter^ (Ar^, Av^) must be inferred from the impact of acceleration on 
the second derivatives of s and a. There is no particular reason to express the final two parameters 
as (Ar^, Avz)', they might, for example, be just as well written as the instantaneous angular acceler¬ 
ation. Regardless, because microlensing events last only a small fraction of an orbital period, clearly 


^Note that as discussed in Appendix A of ISkowron et all 1 2011 1. the microlensing « direction points toward the 
observer, i.e. +z is a blueshift, which is opposite the convention for RV. 
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any such measurements must be substantially less precise than the other parameters. Nevertheless, 
all 13 parameters have been well measured in at least 4 cases (OGLE-2005-BLG-018, OGLE- 


2009-BLG-020.MOA-2011-BLG-09 0. OGLE-2nil-BLG-n417: IShin et al.ll2nill : ISkowron et al.ll2nil 
Shin et ah 2012 : Gould et ah 2013 b 


We now explain how these light curve parameters can be transformed into physical properties. 
With the exception of and the sign of Ar^, all 13 binary parameters can be recovered with 
various combinations of the known parameters. 

As in RV, 6 is automatically measured. In addition, the measurement of p enables a determi¬ 
nation of 0E- This is because the fit to the light curve yields the source flux and so, if there are 
measurements in two bands, the position of the source on an instrumental color-magnitu de diagram 


and h ence the dereddened flux F and surface brightness S', and so finally 0* = y^F/vrS (jYoo et al 


20041 1 . Then 0 e = which allows the transformation of three others into more familiar form 

A6» = S0E, (2) 

where A0 is the instantaneous angular separation between the two components. 


f^geo — 

EE 


The combination of tte and 0 e then add three more binary parameters. First, of course, tte 
immediately yields the direction of /igeo- Then, because 


el = KMlT.eh 


K = 


4G 

AUc2 


~ 8.1 


mas 

Mq’ 


( 3 ) 


(|Einsteinlll936l l and because is usually known quite well, we obtain the mass and distance of the 
system 

M = ^; Di= . (4) 

KTTe C/eTTe + VTs 

The measurement oi Dl enables us to transform angular measurements into physical measure¬ 
ments, i.e., Arx = D^AO, and = DiO. Hence, we now have 9 parameters, 

(mi, m 2 , r, vj_, Ar_L) where r is the 3-space position of the lens system at to- Finally, via the other 
measurement parameters, these yield Av_l = DlsOe'J and so 13 parameters (mi, m 2 , Dl, r_L, vj_, 
Ar_L, Av_l, TAr^, Av^). 

Table [3] summarizes the microlensing parameters, and Table HI gives an overview of h ow those 
parameters translate into the parameters of a binary. Appendix B of ISkowron et al.l (120111 ) provides 
additional details on the transformation between microlensing parameters and the parameters of a 
Keplerian orbit. 


2.3. RV vs. Microlensing: Points of Comparison 


To understand the conditions under which binary microlensing observations can be tested by 
RV, we now focus on 10 parameters, namely the six phase-space coordinates of internal motion 
(Ar,Av), the two masses ((mi,m 2 ), or equivalently {M,q)), the system distance [D^ or r^) and 
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the direction of transverse motion (fi/fj,). We therefore ignore the system radial velocity (vz), which 
can be measured very well by RV but not at all by microlensing, the magnitude of the proper motion 
{fi), for which the reverse holds true, as well as the system angular coordinates (0), which were 
included only for completeness. 

RV observations measure 5 combinations of these quantities, namely four internal phase-space 
coordinates and the mass function, which is a combination of mi, m 2 , and i. Here, we ignore for 
the moment the fact that spectroscopic measurements also return an estimate of the primary mass 
mi- 

We begin the analysis of microlensing by examining the “typical good case”, in which (q, s, 6-^, 
tte, 7 ) are measured, while (Ar^, Avz) are not measured. These eight measured quantities are all 
combinations of the ten parameters under consideration. That is (g,7rE,0E) Dl, fi/fj,), 

and (Ar_L, Av_l) -f-?- Dl6-^{s, 57 ). Hence, with a total of 8 microlensing plus 5 RV = 13 constraints 
on a total of 10 parameters, there are nominally three independent points of comparison. 


However, in the case of OGLE-2009-BLG-020, the situation is not quite so favorable. There is 
a well-known degeneracy between 7 _l and tt e i , the component of tte parallel to the projection of 
Earth’s acceleration on the plane of the sky (IBatista et al.ll2nilh. a nd this degeneracy is present in 
the solution of OGLE-2009-BLG-020 as well ( Skowron et al. 2011 ). Hence, in fact, there are only 
two independent points of comparison in the present case. 


In addition, if we consider that RV is returning six of the 10 parameters (i.e., including the 
spectroscopic determination of the mass) then the parameter counting yields 6 - 1-8 — 1 — 10 = 3 
independent points of comparison. In the present work, we will not consider the mass test due to 
the low signal-to-noise of the spectra. Nevertheless, it is clear that with by combining microlensing 
and RV observations, the problem is over-constrained, allowing a direct test of the microlensing 
orbit prediction. 


3. Observations 


3.1. Microlensing Data 


The microlensing data on OGLE-2009-BLG-020 are described in detail in Section 2 of ISkowron et al 


(j201ll ) and shown in their Figure 1. In brief, the event was monitored by the Optical Gravitational 
Leasing Experiment (OGLE) and Microlensing Observations in Astrophysics (MOA) survey groups. 
In addition, it was observed by eight followup telescopes: Bronberg Observatory 36cm, Campo 
Gatino Austral Observatory 40cm, CTIO/SMARTS 1.3m, Farm Cove Observatory 36cm, Faulkes 
Telescope North 2.0m, Faulkes Telescope South 2.0m. Kumeu Observ atory 36cm, and University 
of Tasmania 1.0m. As described in Section 2.2 of ISkowron et al.l (|201ll ). the data have had outliers 
removed, been binned, and had the error bars rescaled so that y^/dof ~ 1 for each data set. 
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3.2. Keck/HIRES 


Between Mar 2011 and Oct 2011 OG LE-2009-BLG-02 0 was observed six times by the High- 
Resolution Echelle Spectrometer (HIRES; IVoet et alJll994l ) on the Keck I telescope with an expo¬ 
sure time of 700s. One additional observation was taken in Aug 2013 with an exposure time of 900s. 
All the spectra were taken without the iodine cell and R ~ 55,000. A summary of observations is 
given in Table El 


These data were reduced using the standard Galifornia Planet Search (CPS) pipeline (jHoward et al 


20101 1. 


3.3. Magellan/MIKE 


Seven of the observations were taken with the r ed camera of the Magellan Inamori Kyocera 
Echelle spectrometer (MIKE; iBernstein et al.ll2003l ) on the 6.5m Magellan/Clay telescope. The 
observations are summarized in TableEl Most of the observations were taken with the 0.7^' slit with 
a 900 second exposure time. The exceptions were observations on 24 May 2014 and 5 September 
2014, which were taken in poor seeing conditions {2" and > 3'^ respectively). 


These spectra were reduced using the CarPy MIKE pipeline (jKelson et al.ll2000l : lKelsonll2003l l. 
Each spectrum was reduced individually with the exception of the observations of 5 September 
2014. In that case, we use the MIKE pipeline to stack the four spectra to increase the signal-to- 
noise ratio (S/N). The wavelength calibration was done relative to ThAr lamp observations taken 
before and/or after the science exposures. We removed the blaze by fitting the continuum of each 
order with a (Gaussian)-|-(Constant)-|-(Slope), i.e.. 


f{x) = Pi exp 


-{x-p2f 

2P3 


+ PA+ PbX. 


( 5 ) 


This results in a flatter continuum than the standard GarPy reduction. However, the differences in 
the resultant RVs are minimal. We hnd that our method results in slightly smaller RV uncertainties 
~ 10-20%, but the difference in the measured RVs between the two methods is (< 0.1 km s“^) is 
an order of magnitude smaller than the values of the uncertainties (~ 0.5-1 km s“^). 


4. Radial Velocities 


4.1. RV template 

To determine the radial velocities of the lens primary, we cross-correlate the spectra of OGLE- 
2009-BLG-020 with a template from the CPS database. Because the default CPS pipeline only 


















Table 5. RVs of OGLE-2009-BLG-020 


BJD' 

RV 

(km s“^) 

<7rv 

Observatory 

Instr 

Date'^ 

Slit 

Exposure 
Time (s) 

Resolution 

S/N @ 

6000A 

5635.1293 

-36.23 

0.88 

Keck I 

HIRES 

14 Mar 

2011 

0.861” 

700 

55,000 

10 

5668.0545 

-40.37 

0.90 

Keck I 

HIRES 

16 Apr 

2011 

0.861” 

700 

55,000 

12 

5708.0989 

-42.77 

1.35 

Keck I 

HIRES 

26 May 

2011 

0.861” 

700 

55,000 

8 

5723.9147 

-44.12 

0.93 

Keck I 

HIRES 

11 Jun 

2011 

0.861” 

700 

55,000 

10 

5797.7612 

-43.82 

0.74 

Keck I 

HIRES 

24 Aug 

2011 

0.861” 

700 

55,000 

13 

5843.7469 

-37.06 

0.27 

Keck I 

HIRES 

9 Oct 

2011 

0.861” 

700 

55,000 

8 

6530.7713 

-42.96 

1.63 

Keck I 

HIRES 

26 Aug 

2011 

0.861” 

900 

55,000 

12 

6586.4889 

-44.38 

0.85 

Magellan/Glay 

MIKE 

20 Oct 

2013 

0.7” 

900 

31,000 

10 

6586.5058 

-44.30 

0.52 

Magellan/Glay 

MIKE 

20 Oct 

2013 

0.7” 

900 

31,000 

12 

6721.8966 

-33.81 

0.99 

Magellan/Glay 

MIKE 

4 Mar 

2014 

0.7” 

900 

31,000 

15 

6722.8966 

-33.23 

0.56 

Magellan/Glay 

MIKE 

5 Mar 

2014 

0.7” 

900 

31,000 

17 

6801.8221 

-42.13 

0.98 

Magellan/Glay 

MIKE 

24 May 

2014 

1 . 0 ” 

900 

22,000 

7 

6906.5487^ 

-42.83 

0.67 

Magellan/Glay 

MIKE 

5 Sep 

2014 

0.7” 

6900 

31,000 

7 

6908.6027 

-43.25 

0.44 

Magellan/Glay 

MIKE 

7 Sep 

2014 

0.7” 

1800 

31,000 

25 


®'Start of night 

^This observation is the sum of 4 exposures (3 x 1800s and 1 x 1500s). The quoted signal-to-noise 
is for the summed spectrum. 
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produces a flattened and stitched spectrum of HIP63762 covering the wavelength range 5010A- 
6309A, we use only the orders of the HIRES spectrum of OGLE-2009-BLG-020 that overlap with 
that region (i.e. orders 0-14). We combine the indivi dual GGFs for each order to create a maximum 
likelihood function (ML) following the prescription in lZuckerl (j2003l L This produces a double-peaked 
ML in which the taller peak corresponds to the (brighter) lens and the smaller peak to the source 
(see Figured]). 


We selected the template by cross-correlating the highest S/N Keck spectrum (from BJD' = 
5797.7612) with spectra of all stars in the GPS database and measuring the heights of both the 
source and lens peaks. We find that the lens is broadly consistent with being an early K-dwarf and 
select our template from amongst the 10 best matches with 5.5 < My < 8.5 and 0.85 < {B — V) < 
1.35 (cuts that encompass most stars with likelihood peaks within 95% of the maximum value). 
Note that our template choice, HIP63762, is also broadly consistent with the source spectrum, 
obviating the need for a two template fit. Because the source is a clump giant and the GPS sample 
specifically selects against giants, it is unsurprising that the best available template for the source 
would be a dwarf of similar color. 


4.2. Keck Velocities 

To determine the radial velocity of the lens and source for each epoch, we calculate the ML of 
OGLE-2009-BLG-020 as compared to HIP63762 in the manner described above. We then fit the 
two tallest peaks of the ML with 3, 4, and 5 parameter Gaussians (see Equation [5|) , and take the 
result with the best x^- The measured lag of each peak is then the value of p 2 for the best fit. 

Because the source star has a constant velocity (confirmed by applying a barycentric correction 
to the measured lag), we can use it as a wavelength reference in a manner analogous to an iodine 
cell. Both the source and the lens spectra encounter the same systematic effects as their light passes 
through the telescope and instrument optics. Hence, we take the measured RV of the lens to be 

KVlens ^^Ssource ^^Slensi (®) 

where lagsource and lagiens are the lags of the source and lens, respectively, measured relative to 
HIP63762. This automatically takes into account (and removes) the barycentric velocity and 
other systematics induced by the instrument. The uncertainty of the RV is taken to be the 
standard deviation of the RVs measured for each individual order using the same method (i.e., 
RViens,i = lagsouj.ce,i “^agig^g j). The Uncertainties in the RVs are dominated by uncertainties in the 
measurement of the source peak. These uncertainties are larger than they would be if we simply 
made a barycentric correction to the measured velocity of the lens peak. However, the RVs using 
our method are sufficiently precise for characterizing the system and likely to be more accurate 
because our method automatically accounts for systematic effects. The final velocities are given in 
Table [5] where BJD' = BJDtbd — 2450000. 





- 14 - 


4.3. Magellan Velocities 


For the Magellan/MIKE spectra, we used only the orders from the red camera that overlapped 
with the iodine region, i.e. orders 57-68 (12 total). We also cross-correlate those orders against the 
flattened Keck/HIRES spectrum of HIP63762. We create the ML and extract the lag of tallest peak 
in the same way as for the Keck data. However, as can be seen from EigurelH the ML (and CCEs) for 
the Magellan/MIKE spectra are lower resolution and lower S/N than the Keck/HIRES spectra, and 
so we are not able to reliably extract the lag of the source using this met hod. Instead, we calculate 
the barycentric correction explicitly for each observation using BARYCORR (IWright fc Eastmanll2ni4l ) 
and apply it to the measured lens RV. This leaves a velocity offset between the Magellan RVs and the 
Keck RVs equivalent to the difference in radial velocity between the source star and the template. 
We can place the Magellan lags on the same velocity scale as the Keck RVs by calculating this 
difference, i.e. the weighted mean lag of the source star in the Keck data (< laggom-ceK >) after 
the barycentric velocity has been removed. Hence, the Magellan velocities on the Keck system are 
given by 

RMlenS < laggQm-j-g > (lagjg^g ]y[g^g BC). (7) 


To compute the uncertainty in the radial velocities from the Magellan data, we first compute 
the uncertainty in RViens hy computing the standard deviation of RViens measured for each order 
individually, as we did for the Keck data. Then, we add this in quadrature to the standard deviation 
of < laggom-ce Keck > measured from the Keck MLs ('7<iagsource> = 0.29kms“^). Einal values for 
the radial velocities are given in Table [5j 


4.4. RV Orbit 


We use EXOFAST ([Eastman et al.ll2013l i to find the best-fit orbit to the radial velocity data of 
OGLE-2009-BLG-020. This package finds preliminary solutions using a Lomb-Scargle periodogram, 
refines them using an Amoeba minimization, and determines the uncertainties using a Markov Ghain 
Monte Carlo. We began by finding a preliminary solution for the period using a Lomb-Scargle 
periodogram and then seed EXOFAST with this value as a prior. We allow for free eccentricity but 
do not allow for a slope, such as might be caused by a third body in the system. EXOFAST uses 
BJD as the time standard. 


These fits clearly indicate that our uncertainties for the radial velocities are over-estimated, 
and rescales them by 0.4255. This £t shows that the lens has a period of ~ 276 days and an 
eccentricity of 0.341. The red line in Eigure [2] shows the best-fit RV curve to the data. The full 
RV solution (parameters and their uncertainties) is given in Table El The posteriors are shown in 
the lower-left panels of Eigure El Note that EXOFAST provides the argument of periastron of the 
primary cj*. We report the argument of periastron of the secondary Wpeo = w* — 180deg. 
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Keck Magellan 




-300 -200 -100 0 100 -200 -100 0 100 200 

Lag (km/s) Lag (km/s) 


Fig. 1.— Typical maximum likelihood functions for the spectrum of OGLE-2009-BLG-020 from 
Keck/HIRES (left) and Magellan/MIKE (right). The taller peak corresponds to the lens primary 
and the fainter peak to the source star. The velocity scale is arbitrary. 


Table 6. RV Parameters: Median values and 68% conhdence interval for OGLE-2009-BLG-020 


Parameter 

Units 

RV-only 

RV Orbit Parameters: 


^peri. 

Time of periastron (BJDtdb) • • • 

6142.71^1 

p . 

Period (days). 

276.37;j]|^ 

e. 

Eccentricity. 

n QQc:+0.074 

U.OOO_Q Q5g 

i^peri. 

Argument of periastron (degrees) 

156.8 ±3.4 

M . 

Mass Eunction. 

6.20l^;^^ X 10-3 

Other Parameters: 


a . 

Semi-major axis (AU). 

0.8467 ± 0.0076 

K . 

RV semi-amplitude (km/s). 

D.O^U_q 45Q 
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RVs of OGLE-2009-BLG-020 



BJD’ 


Fig. 2.— Measured radial velocities of OGLE-2009-BLG-020L. The red line shows the best-fit orbit 
to the RV data alone, and the black line shows the best-fit orbit from the joint MCMG including 
both the RV and microlensing data. The gray (and pink) liness show joint (and RV-only) hts that 
are 1, 2, and S-u from the best ht (reflected by the shading). 
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5. Comparison of Independent Fits 


Figure [3] compares the RV constraints on the orbit of OGLE-2009-BLG-020L (Section 14.41) 
derived from EXOFAST to the independent constraints on the orbit from the microlensing light curve 
(jSkowron et al .1120 111 ). The microlensing constraints include the weig hting for the Jacobian , the 
Galactic model, and lens flux as described in Sections 3.4 and 4.2.1 of ISkowron et al.l (|2nill ). For 
the purposes of this comparison, we extrapolate tperi for the RV fit backwards to the time of the 
microlensing observations. In addition, note that microlensing uses HJD as the time standard 
rather than BJD, but this difference is many orders of magnitude smaller than the uncertainties 
in the measured parameters. 


The constraints on the RV paramete rs from the mi c rolens ing light curve are derived from 
the MGMG fits to the microlensing data in ISkowron et al.l (j201ll ). Because of the ± degeneracy in 
(s^, 72 ) measured from microlensing, there is a perfect degeneracy between Wpeu and oopevi — 180deg. 
Since both are equally valid, we plot both solutions in Figure [3] leading to periodic behavior in cuperi- 


Figure [3] clearly shows that the constraints on the orbit of OGLE-2009-BLG-020L from the 
radial velocity data are consistent with the observed properties of the orbit measured from the mi¬ 
crolensing light curve. This is the first confirmation of orbital motion of a 2-body system measured 
from a microlensing light curve. 


6. Joint Fit 


We also perform a joint MGMG fit to the RV and microlensing data to determine the best 
constraints on the physical properties of the OGLE-2009-BLG-02 0L system. The joint MGMG is 
performed in the same parameter space as in ISkowron et al.l (j2nill ). i.e. using to, uqw, tes, i*, '^e,e, 
tte.N) 7 ||, 7 ±) Sz, logq, logic, and 7 ^ as the MGMG parameters. These parameters are similar to 
those in Table O with a few substitutions. In place of tE and p we have teff = 'Uq^e and = ptE- 
In addition, we use the caustic width w instea d of s and uow instead of uq and step in the log of q 
and w. As described in ISkowron et al.l ((20111), this parameterization results in faster convergence 
of the MGMG. 

In addition, we include the angular size of the source, 0*, as a chain parameter. A lthough 
this is an observable quantity (0*^o = 4.45/ras, see Section 4.1.1 of ISkowron et al.ll 201 lh . it has 
some uncertainty (7%) for which we want to allow in the Markov chain. To do this, we allow this 
parameter to float, but we apply a penalty for values that deviate from the observed value, i.e. 


^penalty 


7*,0 


0 . 


*,0 


0.07 


( 8 ) 


Finally, we change the sign of t he RV data to acco unt for the sign difference in the RV and 
microlensing coordinate systems (jSkowron et al.l 120111 1. 
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Fig. 3.— A comparison of the orbit constrai nts for OGLE-2009-B LG-020L from independent fits 
to the microlensing light curve (upper right; Skowron et al. IS) and the radial velocities (Ic 


lower 


left; EXOFAST, Section [43]) . In the upper right panels, the black contours show the RV constraints 
overplotted on the microlensing constraints. The center panels on the diagonal compare the pos¬ 
teriors (gray: microlensing, black: radial velocity ). For the mic rolensing constr aints, the colors 
represent the weights of each link as described in Skowron et al. (Section 3.4 of boill ). For the 
RVs, the colors reflect the likelihood as determined from the of each fit. 
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For each link, the MCMC parameters are converted to binary orbit parameters, which are used 
to generate an RV curve. The degeneracy in ±( 52 , 7 ^) leads to a degeneracy in the sign of the RV 
curve. In addition, the absolute RV offset, 7 , is unmeasured. In order to determine the appropriate 
sign of the RV curve for each MCMC link, we generate the corresponding RV model and fit it to 
the RV data with both signs and optimize for the best value of 7 in each case. We take the better 
of the two hts as the “correct” sign, which determines the sign of {sz,'yz) and also the values of 
and Wperi (see Section 12.2j) . 


The results of t he MCMC are shown in the lower left panels of Figure H] in comparison to 
the constraints from ISkowron et al.l ([20111). This clearly shows that including the RV data vastly 
improves the constraints on the orbital solution. To determine the final parameters of the system 
an d their uncerta.inties , we weight the MCMC chain from the joint fit by the Jacobian (Appendix 
B ISkowron et al.ll 2 nilh to account for the transformation from the MCMC parameters to physical 
parameters. The hnal values for the binary orbit are given in Table [71 


Note that ISkowron et al.l (|201lh require that the parameters of the lens star (flux, mass, and 
distance) are consistent with theoretical isochrones. This sets the upper and lower boundaries in 
Mtot for the microlensing-only chain. We do not include this weighting in our joint hts, which is 
why the posteriors extend to regions ’’excluded” by the microlensing MCMC. 


7 . Conclusions 


We have performed the hrst test of a microlensing detection of lens orbital motion by direct 
comparison of the microlensing orbit constraints to the measured orbital parameters from radial 
velocity observations of the lens system. Although the source and lens are not resolved, we show 
that the “contamination” of the lens spectrum by the source star is actually helpful. The fact that 
source is moving at constant radial velocity allows it to serve as a wavelength reference for our 
high-resolution spectra of the lens. We hnd that the orbit of OGLE-2009-BLG-020L as determined 
from radial velocity is fully consistent with the constraints from the microlensing light curve, which 
makes it the hrst conhrmation of a microlensing measurement of orbital motion. 

Demonstrating that the parameters of the microlensing solution are consistent with RV follow¬ 
up is a very strong conhrmation of the method for including orbital motion in microlensing analysis. 
This test is completely independent of the microlensing data and is stronger than many previous 
tests of microlensing results because it constrains more parameters. Hence, we can now view the 
entire microlensing orbital motion sample and the parameters we have derived with more conhdence, 
including in the case of planetary microlensing events for which such followup is not possible. 


In the future, a stronger test should be possible for the microlens OGLE-2011-BLG-0417 
([Gould et al.l l2013l : IShin et al.l l2012l l. While the lens in this event is fainter, the microlensing 
constraints on the orbit are much better. In particul ar the form of the RV curve is predicted from 
the microlensing orbit measurement (see Figure 1 of iGould et al.ll2013l ). 
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Fig. 4.— Constraints on the orbit of OGLE-20 09-BLG-020L. The upper right-hand panels show 
the constraints from the ISkowron et al.l ([20111) MCMC fit to the microlensing data. They are 
reproduced as the gray shaded regions in the lower left panels. The colored points in the lower 
left panels show the results of a joint MCMC to both the RV and microlensing data. The best 
joint fit is indicated by the black circle in the upper panels. The panels on the diagonal show the 
marginalized distributions of each parameter for the microlensing-only fit (gray) and the joint fit 
(black). 
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Table 7. 


Joint Fit: Orbit Parameters 


Parameter 

Units 

Value 

Uncertainty 

Eccentricity 


0.265 

±0.025 

a 

AU 

0.865 

±0.024 

tperi 

days 

4758.692 

±2.685 

^node 

deg 

-7.767 

±1.260 

inclination 

deg 

129.424 

±1.238 

^peri 

deg 

151.600 

±3.325 

Mtot 

Mq 

1.132 

±0.097 

Dl 

kpc 

0.747 

±0.020 

Q 


0.275 

±0.003 

Derived: 




Period 

days 

276.555 

±0.300 

mi 

Mq 

0.244 

±0.021 

m 2 

Mq 

0.888 

±0.076 







